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ABSTRACT 

We present an analysis of a mosaic of ROSAT PSPC and Einstein Observatory IPC 
X-ray observations of 14 clusters and 2 groups of galaxies, enclosed in a sky area 
15° x 20° centred on A3558 within the Shapley Supercluster region. 

From the mass of each cluster, extrapolated to a density contrast of 500, we 
define 4 large structures: (1) the core of the Shapley Supercluster with radius 13 /i^q 1 
Mpc; (2) the core plus A1736; (3) the core, A1736 and the western extension (A3528- 
A3530-A3532) and A3571; and (4), adding the northern cluster pair A1631-A1644 to 
structure (3), the Supercluster as a whole enclosed within a radius ~ 90/i^q 1 Mpc. The 
observed total masses range between 3.5-8.5 xlO 15 M0. The mass values derived from 
the observed intracluster gas (assuming a baryon fraction consistent with primordial 
nucleosynthesis) are of the order of few times 10 16 Af Q . Given these estimates, the 
core is a bound structure with a very significant overdensity of at least 5 times the 
critical density, indicating that is approaching maximum expansion before collapsing. 
Structure (2) has an overdensity of 1.7 on a scale of ~ 30h^ Mpc. The core is then 
a 2-3.7 a fluctuation in an initial Gaussian perturbation field, normalized for galaxies 
and clusters. The highest value applies if we extrapolate to a density contrast of 200 
and assume that fib < 0.095. 

The baryon fraction of the core of the Shapley Supercluster, with A1736, is about 
15 per cent over a radius of 28 h^Q Mpc. 

Key words: galaxies: clustering - X-ray: galaxies. 



1 INTRODUCTION 

First noted by Shapley (1930), the Shapley Supercluster is 
a large overdensity of galaxies in the region of Centaurus- 
Hydra (a 2 ooo : 13 h 25 m , <5 2 ooo : -30°, z ~ 0.046) which has 
been widely studied both in the optical (Melnick & Moles 
1987; Vettolani et al 1990; Raychaudhury et al 1991; Quin- 
tana et al 1995) and X-ray wavebands (Day et al 1991 based 
on GINGA data; Breen et al 1994 on Einstein Observatory 
IPC data). 

Analysis of the structure of the supercluster, performed 
by Zucca et al (1993; SC 26 in their catalogue) and Einasto 
et al (1994; SC 80), suggest that it is a concentration 
with more than 20 clusters of galaxies within a volume of 
60 x 100 x 200 Mpc 3 (~ 50 times the average density of 
ACO clusters at the same \b\ ~ 30°, Vettolani et al 1990; we 
assume Ho = 50 kms -1 Mpc _1 and qo — 0.5), distributed 
in 4 different regions: (i) the core around A3558 and A3562; 
(ii) an eastern part centred around A3571; (iii) a western re- 
gion with A3532 and (iv) an elongation to the north (A1736, 
A1644). The unusual nature of this region is also borne out 
through X-ray observations which show these six clusters 



to be amongst the 46 X-ray brightest at high galactic lati- 
tude (|6| > 20°, Edge et al 1990). Furthermore, the Shapley 
concentration probably contributes about 10-20 per cent to 
the optical dipole observed in the motion of the Local Group 
with respect to the Cosmic Microwave Background (Lynden- 
Bell et al 1988, Raychaudhury 1989, Scaramella et al 1989). 

The Shapley Supercluster is thus an extremly dense re- 
gion where peculiar and global cluster-supercluster char- 
acteristics can be studied. Following an initial analysis by 
Fabian (1991), which showed that its central region is the 
largest mass overdensity yet discovered on so a large scale, 
we attempt here to estimate the matter (gas, stellar, gravita- 
tional) distributions, the total baryon fraction, and whether 
the region is bound. To this aim, we collate all the available 
ROSAT PSPC images of the clusters present of the Shapley 
area, completing the resulting mosaic with Einstein Obser- 
vatory IPC images of 3 other clusters. Then, their surface 
brightness profiles are deprojected to estimate the gas and 
gravitational masses, as described in Section 2. Combining 
these results with constraints on the mass in galaxies, we 
obtain limits on the baryon fraction of the Shapley Super- 
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Figure 1. The spatial distribution of the 47 PSPC circular fields, 
each with a diameter 2°, of the 3 IPC 75'-square field and the 
positions of the ACO clusters are shown enclosed within a radius 
of 15° from 13 h 28 m and -25° 15'. The dashed line indicates the 
region shown in Fig. 2. 

cluster in Section 3. In Section 4 we discuss the matter (gas 
and dark) distribution in the various Shapley structures, and 
the implications of our results for cosmology and large scale 
structure. We summarize our conclusions in Section 5. 



2 X-RAY DATA AND ANALYSIS 

The surface brightness profiles from R OS AT Position Sensi- 
tive Proportional Counter (PSPC) images of 11 clusters and 
2 groups, and Einstein Observatory Imaging Proportional 
Counter (IPC) data for 3 other clusters (A1631, A1644, 
A1736), are analysed using the image deprojection proce- 
dure, as described below. 

The spatial position of the PSPC and IPC fields and 
the clusters of the ACO catalogue suspected as members of 
the Shapley Supercluster are shown in Fig. 1; a mosaic of 
the exposure-corrected images is shown in Fig. 2. 

In Table 1 the ACO coordinates and details of the X- 
ray exposure for each cluster in the sample are summarized. 
The optical data on redshift z and velocity dispersion a, 
are taken from the literature; the X-ray data includes the 
instrument, the exposure time, and the spatial average in- 
tracluster medium (ICM) temperature, Tx, (available from 
previous work). In the last column, we quote the Galactic 
absorption taken from the 21-cm determinations by Stark 
et al (1992). 

2.1 Deprojection analysis 

We use PSPC images in the 0.4-2 keV band, which maxi- 
mizes the signal-to-noise ratio due to reduced background. 



The images are subdivided into four energy bands (42-51, 
52-90, 91-131, 132-201 PI channels) and corrected by the 
corresponding exposure maps using an implementation of 
the procedure written by Snowden and collaborators (1994) 
in the Interactive Data Language. The 0.8-3.5 keV IPC band 
images (corresponding to the range 5-10 PI channels) are 
also flat-fielded. 

The cluster surface brightness profiles are extracted us- 
ing the XIMAGE software (vers. 2.53), and then analysed 
using the deprojection technique pioneered by Fabian et al 
(1981). As a more detailed description is given by White 
et al (1996), we only summarize the assumptions and the 
limitations of this analysis method here. Assuming spher- 
ical geometry for the cluster and hydrostatic equilibrium, 
the deprojection technique enables the volume count emis- 
sivity of the hot ICM to be determined as a function of 
radius from the surface brightness profile. Including the ef- 
fects of absorption by the intervening material (hydrogen 
column density A^h) and detector characteristics, a choice 
of the gravitational potential (functional form, core radius 
r c , velocity dispersion cr dpr ) allows the physical properties of 
the ICM to be determined (assuming the perfect gas law). 
In this analysis, we use a true isothermal sphere to describe 
the shape of the gravitational potential (Binney & Tremaine 
1987; a comparison of the effect of different potential laws 
is described by White & Fabian 1995). 

Given the complexity of the Shapley region, some uncer- 
tainties in the available values of a from optical analysis may 
arise from diffuse clusters overlapping and/or the ill-defined 
boundary of each structure. Because of this, we interpolate 
the Odpr values from Tx where available, using the general 
relationships from White et al (1996) (see notes of Table 1; 
sometimes there can be a significant discrepancy between 
the interpolated er dpr and the optical velocity dispersion, i.e. 
A3530 and A3532). The optical velocity dispersion is used 
when T x is unknown (i.e.for A3556, A3559, A3560, SC1327- 
312 and SC1329-313). The consistency of the optical value is 
then compared against the values inferred by interpolating 
the deprojected Tx and luminosity Lx- If there is a large dis- 
agreement, then <T dpr is altered according to the interpolated 
values (as found for A3559, A3560, SC1329-313). 

In Table 2, the input values of r c , the velocity disper- 
sion adpr and the outer hydrostatic pressure, P ou t, are pre- 
sented. These parameters enable us to produce a flat depro- 
jected temperature profile for each cluster, according to the 
isothermal assumption. The table also shows the results of 
the deprojection analysis. The rebin factor gives the degree 
of grouping of the surface brightness profile to improve the 
counts statistics, and the radius R out is the outer extremity 
of the deprojection where the luminosities and masses (the 
gravitating mass, Md pr , the gas mass, M gas , and the stellar 
mass, M*) are quoted. The values with the subscript "500" 
are discussed below in Sect. 3. 

Note the background contribution is estimated by con- 
sidering the surface brightness value just outside the maxi- 
mum radius of each deprojection. 

3 THE BARYON FRACTION 

Comparison of observations of light-element abundances 
with standard primordial nucleosynthesis theory indicates 
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Figure 2. In the panel below, we present 
a mosaic of the X-ray observations. In 
the panel to the right, we show the obser- 
vations of clusters in the northern exten- 
sion (A1631 and A1644) which is located 
~ 16° away from the core in the NW 
direction. Note the two groups SC1327- 
312 and SC1329-313 (west and cast, re- 
spectively) between A3558 and A3562. A 
close-up image of the core is shown in 
Fig. 7. 
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that the baryon density parameter fib lies between "reason- 
able" values of 0.037 and 0.088 (see Copi et al 1995 for 
a discussion on the meaning of "reasonable" due to the fact 
that the main uncertainties are not described by a Gaus- 
sian). These limits represent fi b ow and fi|j lsh respectively, 



that we consider hereafter when referring to the baryon den- 
sity parameter from primordial nucleosynthesis. It should be 
noted that these changes, from the previously accepted value 
of 0.05 ± 0.01 (Walker et al 1991), arise from revised con- 
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Table 1. The X-ray cluster sample 



name 


inst. 


exp 


<*2000 


^2000 


z 




ref 


T l 


ref 


N H 






sec 


hh mm ss 


dd mm ss 




km s~ 1 


opt 


keV 


X-ray 


10 20 cm- 2 


A1631 


IPC 


5705 


12 52 50 


-15 26 17 


0.0466 


628 [2.7] 


(1) 


2.8 [653] 


(4) 


3.9 


A3528 


PSPC 


15751 


12 54 18 


-29 01 16 


0.0521 


864 [4.8] 


(2) 






6.1 


A3530 


PSPC 


8620 


12 55 37 


-30 21 14 


0.0532 


391 [1.1] 


(2) 


3.2 [ 703] 


(4) 


6.0 


A1644 


IPC 


11096 


12 57 15 


-17 22 13 


0.0475 


933 [5.6] 


(1) 


4.7 [869] 


(5) 


4.3 


A3532 


PSPC 


8620 


12 57 19 


-30 22 13 


0.0537 


594 [2.4] 


(2) 


4.4 [838] 


(6) 


6.2 


A3556 


PSPC 


13765 


13 24 06 


-31 39 38 


0.0481 


554 [2.1] 


(3) 






4.7 


A1736 


IPC 


10660 


13 26 52 


-27 07 33 


0.0446 


835 [4.5] 


(2) 


4.6 [858] 


(5) 


5.1 


A3558 


PSPC 


30213 


13 27 55 


—31 29 32 


0.0475 


986 [6.1] 


(31 


6.2 [1012] 


(7) 


4.4 


A3559 


PSPC 


8127 


13 29 54 


-29 31 29 


0.0461 


279 [0.6] 


(2) 






4.5 


A3560 


PSPC 


14929+ 


13 31 51 


-33 13 25 


0.0462 


403 [1.2] 


(2) 






4.3 


A3562 


PSPC 


20202 


13 33 32 


-31 40 23 


0.0491 


825 [44] 


(2) 


3.8 [772] 


(6) 


4.2 


A3564 


PSPC 


8011 


13 34 22 


-35 13 21 


0.0498 




(2) 






4.5 


A3566 


PSPC 


7357 


13 38 59 


-35 33 13 


0.0497 




(2) 






4.3 


A3571 


PSPC 


6072 


13 47 29 


-32 51 57 


0.0391 


1060 [7.0] 


(2) 


7.6 [1133] 


(6) 


4.0 


SC1327-312 


PSPC 


50415+ 


13 29 45 


-31 36 12 


0.0495 


676 [3.1] 


(3) 






4.3 


SC1329-313 


PSPC 


20202 


13 31 36 


-31 48 45 


0.0469 


1044 [6.8] 


(3) 






4.2 



(1) Zabludoff et al 1990; (2) Quintana et al 1995; (3) Bardelli et al 1994 

(4) White et al 1996; (5) David et al 1993; (6) Edge et al 1990; (7) Day ct al 1991 

t the exposure is equal to the sum of the overlapping images 

t in the square brackets and italic font, we report the reference values (cr rc f, T re f) obtained from 

the relationships (cf. White et al 1996): T rof = (0.093 ± O.OeO)^' 1 ' 831 * ' 291 ) and a rci = (3.697 ± 0.547)t4°' 552±0 ' O88) , 
where a (measured in 100 km s -1 ) and Tx (in keV) are the observed values. 



A1631 scq. 1900; 

A3528 scq. wp300093; also known as Klcmola 21 (Klcmola 1969), it is a X-ray double cluster: 

A3528n (12 h 51 m 38.2 s , -28°44'01", T x = 3.8 keV [772 km/s] ref. 4); 

A3528s (12 h 51 m 58.5 s , -28°57'20", T x = 4.0 keV [795 km/s] ref. 4). 
A3530 seq. wp701155n00; its ROSAT field contains also A3532. 
A1644 seq. 7654; it is a cD cluster. 
A3532 scq. wp701155n00; also Klemola 22. 
A3556 scq. wp800375n00 

A1736 scq. 7653; a foreground group of galaxies is also seen optically at ~ 10400 km/s (Dressier & Schectman 1988) 
A3558 scq. wp800076; Shaplcy8, a cD cluster (cf. also PSPC analysis by Bardelli et al 1995) 
A3559 scq. rp800285n00. 

A3560 seq. wp800381+wp800381a01+wp800381a02. 
A3562 seq. rp800237n00. 

A3564 scq. rp800288n00+a01; no significant X-ray emission; not considered in the sample analysed. 
A3566 scq. rp800286; no X-ray significant emission; not considered in the sample analysed. 
A3571 scq. rp800287; it is a cD cluster. 

SC1327-312 seq. wp800076+rp800237n00; singled out optically in (3) as SC1329-314 B and in the X-ray band 
by Breen et al. (1994). 

SC1329-313 scq. rp800237n00; singled out optically in (3) as SC1329-314 A and in the X-ray band 
by Breen et al. (1994). 



straints on the D + 3 He (for the lower limit on Ob) and 7 Li 
abundance (for the upper limit). 

Thus, if regions that collapse to form rich clusters in 
an Einstein-de Sitter Universe retain the same value of Ob 
as the rest of the Universe, then only a few per cent of 
cluster masses can be due to baryons (mostly gas in the ICM, 
but also stars in galaxies). However, observations of clusters 
reveal a baryon fraction, /b (the sum of the gas fraction / gas 
in the ICM and the stellar fraction /* contributed by the 
cluster galaxies), of 10-30 per cent. 

Historically, X-ray observations have always shown a 
relatively high baryon fraction in clusters (e.g. Stewart et 
al 1984), although White & Frenk (1991) highlighted the 
discrepancy when new tighter and lower constraints from 
nucleosynthesis were published. In the Coma cluster, White 
ct al (1993; cf. also Briel et al 1992) found the ratio between 



the gas plus stellar mass to the total gravitational mass to be 
about 0.226 /ig 15 +0.015. Thus, the disagreement between 
the observed baryon density and that predicted by cosmic 
nucleosynthesis appears to be exceeded by a factor of ~ 3 for 
all reasonable values of the Hubble constant (there is only a 
weak dependence on Ho). 

This problem, sometimes known as the "Baryon Catas- 
trophe", is not unique to the Coma cluster as the discrep- 
ancy appears to be similar, if not worse, in many other clus- 
ters (Henry et al 1993, White & Fabian 1995). The baryon 
content of the Shapley Supercluster has also been estimated 
to be significantly high at 18 per cent, in a 37 Mpc re- 
gion around the core (Fabian 1991). Such an overdensity on 
this large a scale means that the severest problem is then 
the accumulation of baryons within a standard cosmological 
model. 
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Table 2. Dcprojcction analysis results 



iiciiric 


(J A 

u dpr 


' C 


P„„t /10 4 

± out/ 


rebin 


^out 


^X,bol 


JVJ dpr 


-^gas 


+ M t 


^500 


M^rm 
JVJ OUU 


fsoo 




km s ^ 


Mpc 


K cm 


x bin t 


Mpc 


10 44 erg s" 1 


1O 14 M0 


10 13 


Mq 


Mpc 


1O 14 M0 


xl00% 


A1631 


628 


0.2 


0.4 


5 


0.40 


0.1 


0.99 


0.17 + 0.42 


1.23 


3.07 


0.060 


A3528n 


772 


0.3 


4.2 


3 


0.37 


1.1 


1.09 


0.62 


+ * 


1.65 


7.23 


0.079 


A3528s 


795 


0.2 


4.0 


3 


0.37 


1.4 


1.40 


0.62 


+ * 








A3530 


703 


0.4 


1.0 


3 


0.76 


1.0 


2.27 


1.61 - 


r 0.50 


1.40 


4.69 


0.117 


A1644 


856 


0.5 


4.0 


3 


0.61 


2.5 


2.10 


1.78 - 


h 0.32 


1.70 


8.12 


0.120 


A3532 


791 


0.4 


1.5 


3 


0.90 


2.9 


3.45 


3.33 - 


h 0.66 


1.55 


6.38 


0.136 


A3556 


644 


0.6 


0.4 


4 


1.55 


0.7 


4.17 


3.82 


+ ... 


1.29 


3.63 


0.102 


A1736 


835 


0.7 


3.5 


3 


0.63 


1.4 


1.62 


1.45 + 0.27 


1.66 


7.54 


0.129 


A3558 


986 


0.6 


2.8 


3 


1.03 


9.8 


5.57 


6.84 + 1.20 


1.94 


12.18 


0.176 


A3559 


511 


0.7 


0.2 


3 


1.19 


0.3 


1.97 


1.89 


+ ... 


1.00 


1.66 


0.104 


A3560 


738 


0.5 


0.7 


3 


1.34 


2.3 


4.67 


5.36 


+ ... 


1.48 


5.43 


0.135 


A3562 


772 


0.6 


1.0 


3 


1.24 


3.9 


4.50 


5.86 - 


h 0.75 


1.54 


6.12 


0.152 


A3571 


1060 


0.3 


4.0 


3 


0.91 


14.5 


6.22 


5.82 - 


r 0.46 


1.99 


12.79 


0.142 


SC1327-312 


676 


0.5 


2.2 


3 


0.60 


1.1 


1.32 


1.26 


+ ... 


1.36 


4.27 


0.137 


SC1329-313 


511 


0.4 


1.0 


3 


0.45 


0.3 


0.61 


0.46 + 0.28 


1.07 


2.05 


0.146 



t the original binsize is 15" in PSPC images and 16" in IPC; 15" ~ 20 kpc at z ~ 0.05. 

A for this double cluster, we estimate that M* = 0.66 X 10 13 Mq at R = 0.74 Mpc, which is the radius that encloses both 
substructures. 



Possible solutions to the "Baryon Catastrophe" are: 
(a) a non-zero Cosmological Constant A, so that fio.tot = 
+ Oa = 1 (cf. Carroll et al 1992 for a discussion on 
the limit on A), (b) a different estimate of fib from the nu- 
cleosynthesis calculations of the light elements abundances 
(as may appear from the recent conflicting observations of 
deuterium abundances, cf. Hata et al 1996 and Schramm & 
Turner 1996, and references there), (c) inappropriate X-ray 
modelling of clusters which may result in overestimations 
of gas mass and/or underestimations of gravitational mass 
(e.g. Gunn & Thomas 1996). 

In our case, the gas fraction, / gas , is estimated by the 
ratio of gas mass determined in the deprojection analy- 
sis to the gravitational mass. To estimate the stellar frac- 
tion, /*, we adopt the optical luminosity function as ob- 
tained by Raychaudhury et al (1991) for the 9 clusters 
and 1 group in common with our sample. Assuming a 
mass-to-light ratio appropriate for elliptical galaxies of 2.46 
x(Lb/1O 1o L q ) o - 35 /i5o(M/Z/) (White et al 1993), we esti- 
mate the stellar mass using the Schechter luminosity func- 
tion (1976) with a slope of a = —1. As the optical luminosity 
function is quoted at 2 Mpc, M* is interpolated to a value 
at -Rout, assuming a King profile for the mass distribution. 

For those 3 clusters (A3556, A3559, A3560) and the 
SC1327-312 group lacking the information required to de- 
termine their optical luminosity functions, we adopt the me- 
dian statistic as a robust estimator of the stellar fractions 
available. The median value of 2.1 per cent (note, the whole 
range is between 0.7-4.6 per cent) is used to determine M* 
in these four objects. The results on / gas and /*, with the 
total baryon fraction /b = / gas x h^" B + /*, appropriate for 
individual values of .Rout, are shown in Fig. 3. 

In this figure, a significant disagreement between the 
maximum value for fib obtained from primordial nucleosyn- 
thesis analysis (dot-dashed line) and that determined from 
X-ray observations of clusters is evident. This is primarily 
due to /gas, which increases with radius, as the contribution 
of /* is almost radially independent and consistent with all 
nucleosynthesis constraints (the dispersion around the aver- 



age value of 0.023 is 0.012). This emphasizes the problem 
that any reasonable / gas , that is ~ 5 ± 3 times /* in our 
sample, would result in a baryon catastrophe. 

In order to compare the fb values within a consistent re- 
gion for all clusters, we extrapolate the matter distributions 
to a density contrast equal to 500. The density contrast, 5r, 
corresponding to a comoving region of radius R, is defined 
as: 

_ (M dpr - Mp) 
Sr - ' W 

where Mo is the expected mass in the same region in a 
homogeneous universe where the density p is equal to p c fio: 

M = |I# 3 "o ~ 2-9 x W^Ri^Qoh^MQ. (2) 

The overdensity value of 500 is chosen as representative 
of the cluster region where the assumption of hydrostatic 
equilibrium used in deprojection is still valid. In fact, various 
numerical simulations (cf. Evrard et al 1995, Bartelmann 
& Steinmetz 1996) show the gas to be largely thermalized 
inside the relaxed regions characterized by this overdensity, 
with outer envelopes that accret and/or merge on the cluster 
at lower density contrast. 

Thus, to determine the baryon fraction /500, we need to 
extrapolate Md pr , M gas and M* to the radius R500, where 
Sr = 500. The resulting values of R500 have a median value 
of 1.54 Mpc and lie between 1.00 (A3559) and 1.99 (A3571) 
Mpc. The extrapolated Md pr , M500 (cf. 11th column in Ta- 
ble 2), is obtained using an isothermal profile with the <7d pr 
and r c as given in Table 2; the baryonic mass at R500 is cal- 
culated (i) by fitting a corrected King model with (3 equal 
to 2/3 to the deprojected density profile and integrating the 
best-fit to obtain the gas mass, and (ii) extrapolating the 
stellar mass as discussed above. The last column in Table 2 
quotes the calculated values of /500 at the same density con- 
trast for each cluster (cf. also Fig. 4). We obtain a median 
estimate of /500 = 0.112/i5 1,5 + 0.018, for the extrapolated 
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Figure 3. The baryon fraction /b (asterisk), the gas fraction / gaa (triangle) and the stellar fraction /* (square), as reported in Table 2, 
are shown for each cluster in our sample. The horizontal dot-dashed line represents the upper limit of 0.088 on Q^, (Copi et al 1995), 



and the dashed line is the radial / gas dependence reported by White & Fabian (1995): / g 



0.0579 + 0.0556H. The errors on / gaa 



come from Monte-Carlo replications of the deprojection determinations (on M gas ); the errors on /* are propagated from the luminosity 
function. In these figures the north and south component of A3528 are considered as one. The apparent deficit of baryons in A3528 is 
due to the different regions used to determine / gas and /*. Values for / gas in both of the northern and southern components can be 
determined individually, while /* is calculated from the optical luminosity using a radius which encompasses both components (see note 
* in Table 2). 



gas and stellar fraction, respectively, with a range extending 
from 0.060 (A1631) to 0.176 (A3558). 



3.1 The Uncertainties in the Baryon Fraction 

The uncertainties inherent in the single-phase deprojec- 
tion analysis are discussed by White & Fabian (1995). We 
note that the method gives direct and tight constraints on 
the gas mass through the comparison with the observed 
emission-weighted temperature. Also, any multiphase anal- 
ysis of the ICM in a cooling flow will not lower the estimated 
baryon fraction at radii larger than the cluster core (Gunn & 
Thomas 1996). Therefore, the main uncertainties arise from 
the determination of the total gravitational mass. 

Our estimates of the total mass as extrapolated val- 
ues at density contrast of 500, M 5 oo, were obtained using 
optical velocity dispersions (and X-ray temperatures) mea- 
sured within iisoo and adopting the isothermal assumption 



valid on the so-defined regions. However, given the complex- 
ity of the Shapley Supercluster because of the great number 
of extended sources present in the field, some uncertainties 
could arise on the total amount of mass present in each clus- 
ter. To limit this uncertainty, we check the cluster masses 
against optical constraints. Quintana et al (1995) have ap- 
plied the Virial Theorem to 8 of the 12 clusters and 1 of 
the 2 groups in our sample (cf. their Table 7 for estimates 
on M g from galaxies within 2 Mpc, which is an acceptable 
limit for the thermalized region as above shown). Using a 
linear dependence of mass with radius, we compare our val- 
ues of M500 with their rescaled M g . Generally, the optical 
value is subject to large errors, both because of the unknown 
shape of the 3-dimensional tensor of the velocity dispersion 
and because of projection effects. In fact (e.g. see David et al 
1995), there may be discrepancies of order a factor 2 between 
the cluster mass determinations from the application of the 
Virial Theorem to optical data and the results from X-ray 
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Figure 4. This plot shows the estimated baryon fraction, /t,, and 
its extrapolated value to 8r = 500, /50O1 f° r the clusters in our 
sample. The dashed lines connect the 2 values for each cluster. 
The vertical dotted line indicates the density contrast equal to 
500. 



analysis. Assigning a gravitational radius equal to R500, wc 
determine a median ratio, with respect to the X-ray value, 
of 1.2 (5th and 95th percentile: 0.7 and 2.3). The extreme 
cases (A3530, A3559, A3560 at the lower end, and A3556 
and A3562 at the upper end) reflect the disagreement be- 
tween the velocity dispersion obtained optically and that 
inferred from Tx. 

Here we note, and then show in Fig. 5, that only by 
a correction factor larger than 1.5 on the gravitating mass, 
can the baryon fraction /500 drop down to the upper limit 
from the nucleosynthesis calculation (1.48=0.130/0.088, for 
h 5 o — 1; for h 5 o > 1 the factor increases till to 3.82 for 
H = 100 km s" 1 Mpc" 1 ). 

In Fig. 5, we consider corrections that could theoreti- 
cally reconcile the estimates of Qh as measured in clusters 
with the value from the primordial nucleosynthesis calcula- 
tions: (a) a contribution to the gravitating mass, although it 
appears unable to make the necessary correction (when it is 
less or equal to a factor 1.5 as discussed above); (b) Oo ~ 0.3, 
which allows agreement within the constraints on Qh , for all 
values of Ho in the range 50-100 km s _1 Mpc -1 ; (c) a com- 
bination of the 2 above corrections, requiring Ho > 78 km 
s _1 Mpc -1 to be consistent with the O{ ) ow . 

However, note that any conclusion on the median val- 
ues observed in our sample become more dramatic when the 
more luminous clusters such as A3558 and A3571 are con- 
sidered. Fig. 4 emphasizes that there is a real dispersion in 
the baryon content of the clusters. 



Figure 5. Constraints on the Hubble constant Ho and the differ- 
ently estimated baryon fraction: (dashed line) the upper and lower 
limit, n{ j ow and n£ igh , on the baryon density parameter from the 
primordial nucleosynthesis calculations of Copi et al 1995; (solid 
line) the median baryonic contribution on the total mass for the 
clusters in our sample, with and without a correction factor of 
1.5 on the gravitating mass and Qo = 1; (dot-dot-dash line) the 
previous case assuming Qo = 0-3. The thick solid line indicates 
the baryon fraction f£ observed in the 3 structures defined in 
Sect. 4. 



4 DISCUSSION 

Using the isothermal gravitational mass profiles, extrapo- 
lated to the density contrast of 500, M500, we attempt to 
find how much of the Shapley Supercluster is bound. Then, 
the baryonic component of the structures are studied and the 
significance of their total mass assessed in terms of Gaussian 
random fluctuations in the universal mass density field. 

In this analysis on the mass component distribution in- 
ferred from X-ray observations of clusters, we use the sample 
defined above and collect all the data relevant to this aim. 
Furthermore, to provide information on the completeness of 
our sample, all the clusters present in the Shapley region, 
as obtained by cross-correlating the lists from Zucca et al 
(1993), Einasto et al (1994) and Quintana et al (1995), are 
binned with respect to their comoving separation R from 
A3558 and shown in Fig. 6. The comoving radius R is de- 
fined as 



R 



£c 2 A35 



+ D 2 r 



2D 



c,A3558 



D Cj i cos 8, 



(3) 



where 9 is the angular apparent separation on the sky be- 
tween A3558 and the cluster i, and D c is the respective 
comoving distance. 

As it is clear from this figure, not all the clusters in the 
Shapley region were detected or observed in X-ray band. 
Within 30 Mpc of A3558, there are 5 clusters undetected: 
3 from PSPC observations (A3554, A3564, A3566) and 2 
(A3552, A3555) from GINGA scan of the region (Day et 
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Figure 6. The histogram of the counts of ACO clusters vs. 
their comoving distance from A3558 is shown, ovcrplotted with 
a thicker histogram of the counts of the X-ray observed clusters 
(also in bold font for emphasis). Within each bin, the clusters are 
labelled in increasing comoving distance, D c , upwards. The red- 
shift measurements are quoted in the present paper and in: (1) 
Quintana ct al 1995; (2) Postman ct al 1992; (3) Abcll et al 1989. 

al 1991). On the other hand, beyond 30 Mpc, 8 out of 10 
clusters were not observed; the other 2, A3572 and A3575, 
were, but remain undetected. 

We also checked for any feature of the X-ray emis- 
sion from these 15 clusters in the ROSAT All-Sky Survey 
Bright Source Catalogue (Voges et al 1996). Within 0.5° 
of the optical centre as reported in ACO catalogue, three 
(A3524, A3557, A3572) have a suspected hard and extended 
counterpart, an other two (A3564 and A3577) have a hard 
but not-extended nearby source. From the respective counts 
rate as quoted in the catalogue, and a conversion factor of 
1.8 x 10 -11 erg cm~ 2 s _1 for 1 count s _1 , we obtain only for 
A3572 a luminosity of about 3.5 xlO 44 erg s _1 (correspond- 
ing to around 2.5 xlO 14 Mq in M gas , when its luminosity 
is scaled to the Coma cluster), whereas the other 4 clus- 
ters have an estimated luminosity of 0.1-0.3 xlO 44 erg s~\ 
equivalent to M gas ~ 0.6-0.9 xlO 14 M©. 

To summarize, we do not expect a significant contri- 
bution from X-ray detectable baryonic matter extra within 
radii less than 30 Mpc from A3558. Some conclusions on re- 
lated topics, such as possible corrections to our estimates of 
the baryonic and total mass on scales larger than 30 Mpc, 
will be discussed in section 4.1 and 4.2, respectively. 

We analyse the distribution of matter that we observe in 
4 hierarchical spherical structures, increasing in size from the 
core (A3558, A3556, A3559, A3560, A3562, SC1327-312 and 
SC1329-313; Fig. 7) to the whole supercluster, of comoving 
radius R centred on A3558 and given by eq. (3) . 

We see that the core, structure (1), has a total mass of 
3.5 x 1O 15 M equal to the sum of M500 for the 5 clusters 



Figure 7. Mosaic of the core of the Shapley Supercluster using 
all the PSPC observations covering the region. 

and 2 groups enclosed within a scale of 13 Mpc. The next 
structure, (2), is defined as the core and A1736, with a total 
mass of 4.3 x 1O 15 M within 28 Mpc. On a scale of 54 Mpc, 
structure (3), which contains the core, A1736, A3571 and 
the western extension (A3528-A3530-A3532), has a total 
mass is 7.4 x 10 5 Mq. Finally, the Supercluster as a whole, 
structure (4), extends to a scale of 88 Mpc, with a total 
mass of 8.5 x 1O 15 M0. Note that these mass estimates are 
highly conservative limits on the total mass in each cluster. 
No account has been made of the mass in the outer region of 
the clusters, where either the isothermal or the hydrostatic 
equilibrium assumption could fail. Extrapolating M dpr to 
8 — 200 (where the cluster reaches virialization) , increases 
the masses by a median factor of 1.49 (in a range 1.40 - 
1.93) with respect to M500, enlarging by the same factor the 
mass estimates for each of the 4 structures. 

A different approach for estimating the expected mass 
in these structures, considered as self-gravitating systems, 
is to apply the Virial Theorem (Heisler et al 1985) on the 
clusters. As shown by Raychaudhury et al (1991), the re- 
gion is retarded from the cosmic expansion through mutual 
attraction. Assuming the extreme condition that each part 
of the Supercluster is virialized, enables us to calculate the 
velocity dispersion as a measure of the gravitational energy; 
we obtain 395, 482, 1212 and 1121 km s~\ for the core, the 
structures (2), (3) and (4), respectively. The harmonic radius 
in each system is also determined, and gives Virial masses 
of 19, 35, 367, 400 xlO 14 M0, respectively. Again it appears 
that the core has enough gravitating mass to account for 
the Virial value, as does structure (2). For the structures 
(3) and (4), the Virial mass is about 5 times larger than the 
respective M grav . 

Adopting the value of ~ 4 x 10 16 M Q as an estimate of 
the mass in the whole supercluster, we find a contribution 
of almost 30 km s _1 to the component towards the Shap- 
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ley Supercluster optical dipole, which is around 520 km s _1 
(Smoot et al 1991). This is less than 6 per cent of this total, 
in a flat Universe with h^o = 1, and requires an increase in 
the cosmological factor hso/Qo' 4 (i.e. a low density universe 
and/or higher Hubble constant) if it is to be more significant 
(cf. Quintana et al 1995). 



4.1 On the distribution of the baryonic matter 

When we consider the sum of the stellar, gas and gravitat- 
ing masses contributed by each cluster (cf. Table 3), then 
for each of the 4 structures we obtain a baryon fraction 
/jf greater than 0.13/i^ 1,5 and, in particular, almost 15 per 
cent within 28 Mpc. Although this result disagrees with the 
predicted range according to the relaxed constraint from 
the standard big-bang nucleosynthesis model adopted here 
(Fig. 5) , it is approaching reasonable agreement with the low 
deterium abundance detection recently reported by Tytler 
et al (1996) of Q b = 0.096±^2o/i,T 2 . To improve on the pre- 
vious estimations of the baryon content for the Shapley core 
(cf. Fabian 1991, Makino & Suto 1993), we place tighter 
limits on the stellar contribution and avoid the use of any 
empirical relationship among the X-ray quantities by mod- 
elling each cluster separately. 

So far we have neglected any contribution from an in- 
trasupercluster medium (ISCM). If we consider the la limit 



cm ~ s 1 ster 1 (from the GINGA Large Area Proportional 
Counter; Day et al 1991), we obtain a lcr upper limit on the 
ISCM mass, Miscm, of 0.5, 1.7, 3.6 and 5.0 xl0 14 M o for 
structures (1), (2), (3) and (4), respectively (using the ther- 
mal bremsstrahlung relation, Miscm oc y- 1 / 4 ^/ 2 j 1 / 2 j anc i 
a King model density distribution with a core radius R c , 
from structure (1), of 13 Mpc and a temperature of 1 keV). 
These estimates are between 1 and 6 per cent of the corre- 
sponding total gravitating mass, and between 11 and 52 per 
cent of the total gas mass observed in clusters. Given this re- 
sult, if the detection of intrasupercluster gas were confirmed, 
the baryon catastrophe in these structures (i.e. on scale en- 
closed between 10 and 90 h^g Mpc) would become more 
severe, increasing the gas component by a factor of 1.1-1.5. 
We note that GINGA cannot place limits on any plausible 
intergalactic medium of temperature much below 1 keV. 

In Fig. 8, we show the dependence with radius, as mea- 
sured in the different structures, of the baryon parameter, 
Ob = Pb/pc = /if ^o[5fl(Af grav ) + 1]. This parameter, on 
scales where the overdensity produced by M grav with respect 
to the background mass M is still important, retains infor- 
mation on both the baryon fraction and the overdensity Sr. 
Hence, the large baryonic overdensity in the Shapley core 
indicates accumulation in a clearly defined potential well, 
whereas a deficiency is evident on a scale larger than 30 
Mpc. It is difficult to accept that there has been accretion 
of baryons from the outer part to the supercluster core (the 
crossing time permits such motion on scales <C 10 Mpc). 
Thus, part of the observed relative deficit of baryons out- 
side 30 Mpc could be due to the incomplete detection in the 
X-ray band of the total gas present. If we consider the me- 
dian values observed in our sample of 6x 10 14 M Q , for M grav , 
and 0.13, for f b , and a further contribution of a factor 1.5 on 
the total mass when extrapolated at S = 200, we are able to 
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Figure 8. The baryonic parameter = pb/pc is shown as func- 
tion of the scale of the 4 different structures, which arc defined 
in Table 3. The two shaded regions indicate the lower and upper 
constraints from the primordial nucleosynthesis as here adopted 
(Copi et al 1995) and from the recent low deuterium abundance 
detection (Tytler et al 1996). H is 50 km s" 1 Mpc" 1 . The er- 
rors, propagated on the uncertainties of the baryon density, have 
dimension comparable to the size of the points. 



increase the fi b values in Fig. 8 to 1.19, 0.23, 0.05 and 0.02 
for structure (1), (2), (3) and (4), respectively. 



4.2 On the cosmological consequences 

In Table 3 the mass estimates in a homogeneous universe 
expected from eq. (2) are compared with the total masses 
extrapolated from the X-ray analysis. We note again that 
our mass estimates using M grav are conservative, since we 
restricted the determinations to Rsoo- A reliable extrapola- 
tion to 7?200 gives an increase of ~ 50 per cent to the cluster 
masses, and thus also to the structures we define. Here, we 
consider another way to constrain the total mass present 
in these structures, assuming that all the observed baryons 
correspond to the fraction expected from the primordial nu- 
cleosynthesis, i.e. 



M" f 

M PN = x n = M grav x x a . 

"b "b 



(4) 



In other words, this is the gravitational mass associated with 
a baryonic mass Mb for an abundance equal to fib in the 
structures defined above (M b is calculated through obser- 
vations of the gas and stellar matter, and the value for O b 
is assumed here to be 0.06, an average value of the 2 limits 
from Copi et al 1995 and within the 2-a range provided by 
Tytler et al 1996). In Table 3 we quote Mpn for each struc- 
ture, assuming Qo = 1. (Again, an extrapolation to 7?200 
increases Mpn by a factor ~ 1.5). 



© 0000 RAS, MNRAS 000, 000-000 



10 



Table 3. Shapley structures observed on a different scale R. For each structure, we quote: in column N c i, the number of clusters analysed, 
the number of clusters without X-ray detection, the number of clusters not observed at all, and, in brackets, the total number of clusters 
in the catalogue (cf also Fig. 6); the radius R where we observe the baryon fraction f£, the ratio between the total baryonic (gas and 
stellar) matter and the total mass M grav = ^ M500 ; the mass Mpn required to observe a baryon fraction equal to the predictions from 
nucleosynthesis; the virial mass M v ; r ; the mass Mo expected in a homogeneous universe; the positive ovcrdensity S R from eq. (1) related 
to the different mass estimated (details in Table 4). The errors on are obtained by propagation of the uncertainties on the stellar 
(errors on the luminosity function) and gas masses (using a 20 per cent uncertainty as upper limit obtained from the extrapolated values 
when compared with the observed ones). Uq is assumed equal to 1. 
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Using Mg rav , M V [ T and Mpn as reasonable mass esti- 
mates, we try to assess the probability of obtaining the den- 
sity fluctuations related to the structures defined above. The 
initial distribution of these fluctuations is assumed to be 
Gaussian in an Einstein-de Sitter universe. Using a spher- 
ical collapse model (cf. Appendix A), the initial comoving 
radius Ro of the sphere that contains at the present time 
the mass observed with overdensity Sr within a comoving 
radius R is given by: Ro = R(l + 5h) 1/3 /(1 + *)■ The 

power 

observed on this scale Ro, written with its dependence on 
fi (Peacock & Dodds 1994, eq. 41), is 



o -0.3 
2 _ "0 



P k W 2 (kR )k 2 dk, 



(5) 



where k is the wavenumber in units of Mpc -1 , W(kRo) is 
a step window function that averages out fluctuations on 
scales smaller than R , and the functional form of P k is the 
CDM-like power spectrum described by Peacock & Dodds 
(1994), normalized to the averaged contribution of different 
classes of galaxy and galaxy clusters. It includes corrections 
for bias, non-linear evolution and redshift-space distortion 
effects and shows good agreement with the datasets and the 
COBE results on scales larger than 20 Mpc. The same 
results in the following discussion can be argued using a 
more physical CDM power spectrum, as approximated in the 
literature (see e.g. Bardeen et al 1986) with a normalization 
of a(16ft 50 1 Mpc)~ 0.6 — 0.7, whereas a forced large-scale 
normalization to the COBE result gives lower values for the 
power by a factor 2, at least on the scales discussed below 
(cf. Peacock & Dodds 1994). 

To compare the CDM-like power prediction with the 
overdensities observed in the Shapley region, we extrapo- 
late from the non-linear overdensity S R to the corresponding 
present linear value S R() (see for details Appendix A). Then, 
the probability that at some point a density fluctuation <5o, 
with respect to the power on the same scale, exceeds the 
observed 5 Rg is given by: 



P(So > S Ro ) = 



exp 



dt, 



(6) 



where t = S /a Ro and v = S Ro /cr R() . 

Using the total mass M grav , the core is the only struc- 
ture with a very significant overdensity of at least 4.05 and 
a probability P of ~ 5 per cent (i.e. 2.0 times a Ro , the 
r.m.s. power value). Using the overdensity upper limit value 
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Figure 9. The total mass as sum of the cluster masses, M grav , 
the expected gravitating mass from the observed baryons, Mpj^, 
and the virial mass, M v ; r , are shown as a function of the radius 
R of the structures which are defined in Table 3 (see text for 
details). The crosses indicate the values extrapolated from the 
galaxy counts (Vettolani et al 1990, Zucca et al 1993), assuming 
that (i) the volumes are centred on A3558, (ii) a representative 
mass for the clusters equal to the median M grav in our sample. 
The dotted line indicates where the density equals the critical 
density for Qo = 1. The errors, with dimension similar to the 
size of the points, come from propagating the respective quanti- 
ties, assuming for the masses the uncertainties estimated in the 
deprojection analysis and for S7 b a dispersion of 0.03, consistent 
with the Copi et al result. 



of 11.47 (from Mpn), structure (1) increases its significance 
to P ~ 0.2 per cent (3.2 times the corresponding r.m.s. 
power). These overdensities suggest that the core region is 
close to the turnaround point, where the perturbed region 
ceases to expand according to the Hubble flow and begins 
its collapse. The perturbation has evolved to its maximum 
expansion in ~ 10 — 13 Gyr (the age of the universe in the 
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present scenario is 13 Gyr) and its collapsing period will last 
for almost 10 Gyr more, under the cosmological assumptions 
discussed above. 

Structure (2) shows an overdensity of 0.66 within a ra- 
dius of 28 Mpc from A3558, based on Mpn. This corresponds 
to a probability P ~ 26 per cent (1.1 x<tr ) and indicates 
structure (2) as a system reaching the turnaround point in 
about 50 Gyr, before collapsing. 

The fio dependence of Mo and Mpn is linear, whereas 
M grav is independent of it (i.e. only Mpn provides a value of 
the overdensity 5r that does not depend upon any assumed 
cosmological model). Thus, in a low density universe (e.g. 
fio ~ 0.3) the overdensities associated with M grav increase 
considerably: in the core and in structure (2) Sr is equal to 
15.8 and 1.3, respectively, when measured with M grav (cf. 
Table 4); structure (3) has an overdensity of 1.6, estimating 
it by M v i r . For fio < 1, the qualitative description discussed 
above on the state of the evolution of the 4 structures is 
still valid, given the collapse condition that S > (fi -1 — 
1)/(1 + z) ~ 2.2 for fio = 0.3 (the corresponding condition 
for fi = 1 is S > 0). 

Neither structure (3) or (4) exhibits a significant den- 
sity peak, when we use fio = 1- Similar results can be ob- 
tained from the catalogue clusters counts (Zucca et al 1993, 
Vettolani et al 1990). They describe the Shapley concentra- 
tion as the richest supercluster known, with an overdensity, 
with respect to the average density of ACO clusters, of ~ 10 
within a region of 1.4xl0 6 /i 50 3 Mpc 3 enclosing 25 clusters, 
to ~ 400 for the 9 clusters contained in the core (as defined 
by Vettolani et al 1990). This overdensity of counts can be 
translated in mass overdensities using the median value for 
Mg rav calculated in our sample, obtaining values in agree- 
ment with our estimates as shown in Fig. 9. 

We now return to the issue of the mass of each cluster 
and of the structures, which are underestimated. Rework- 
ing Table 4 using an extrapolation to R200, we find that the 
core is a 3.7 and 3.5 a fluctuation for fi = 1 and 0.3, respec- 
tively, using Mpn- Only if fib > 0.095 docs it become less 
than a 3 a fluctuation. The 15 undetected clusters (within 
90 Mpc) may provide a further contribution of the order 
of the observed M gra v, increasing by a factor 2 the lower 
limit on the overdensity (using the median value for M500 
of about 6 xlO 14 M0). We emphasize that we are assuming 
the structures to have a spherical shape, which is unlikely for 
large scale collapsing structures (cf. Peebles 1993; see also 
the conclusion on a "cigar-like" shape of the Supercluster by 
Quintana et al 1995 and the recent evidence for a "bridge" in 
the distribution of galaxies between the core and the western 
extension; Bardelli, private communication). 



5 CONCLUSION 

From our analysis of 12 ACO clusters and 2 groups ob- 
served with the ROSAT PSPC and Einstein Observatory 
IPC, we have estimated the baryon and total gravitating 
mass distribution in the Shapley Supercluster. At a density 
contrast, with respect to the expected background matter, 
of 500 we find that the median baryon fraction observed in 
the individual clusters is /500 = 0.11/i^" ' + 0.02. This is 
inconsistent with the standard nucleosynthesis model, un- 
less fio ~ 0.3 — 0.5 or the low deuterium abundance is con- 



Table 4. Significance of the observed positive overdense regions 
(i.e. 5 R > for Q = 1, and Sr > 2.2 for CI = 0.3). The quantities 
here reported are discussed in the text (cf. Appendix A). Starting 
with observed quantities R and Sr , we extrapolate the linear value 
Ro and 8 R(I by cq. (A4-6) for H = 1 and eq. (A12-15) for H < 1. 
The power cfr is obtained in cq. (5). It can be shown that the 
correction to R for f!o 7^ 1 is not significant. Note: the letter a 
means that the corresponding <5r has been calculated using Mpn; 
b, using Mg rav ■ 
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9.7 


(1)6 


13.4 


4.05 


21.9 


1.03 0.53 


13.7 


(2) a 


27.8 


0.66 


31.4 


0.44 0.39 


49.4 




= 0.3 






to = 15.8 Gyr 




(1)6 


13.4 


15.83 


32.8 


1.33 0.45 


13.9 


(l)o 


13.4 


11.47 


29.7 


1.28 0.49 


15.9 



firmed (from Tytler et al 1996, the 1-a upper limit on fib is 
0.120). Cluster to cluster, there are also a significant differ- 
ence in /5oo- In the Shapley Supercluster region, we observe 
a baryon fraction of 15 per cent within 30 Mpc; a deficiency 
beyond this scale is probably due to an undetected warm 
intergalactic medium. 

The gravitating masses M gra v , calculated as the sum of 
the cluster masses within various structures, range from 3.5 
to 8.5 xlO 15 M , which are lower than the minimum neces- 
sary to bind the structures (i.e. Sr > 0, for fi = 1; Sr > 2.2, 
for fi = 0.3) except the core (A3556, A3558, A3559, A3560, 
A3562, SC1327-312, SC1329-313). Applying the Virial theo- 
rem to the structure does not indicate masses equal the crit- 
ical one Mo. Then, on scales of 50 Mpc and larger, i.e. for 
structure (3) (core—west-A1736-A3571) and (4) (the whole 
supercluster) , we compute a mass similar to that estimated 
in a different way by Fabian (1991) and Quintana et al 
(1995), i.e. > 1O 16 M0. Assuming, from primordial nucle- 
osynthesis, a baryon density parameter of 0.06, we expect 
a total mass Mpn of 0.9 xlO 16 M in the core, increasing 
to 1.8 xlO 16 M for the whole Supercluster. These indicate 
that either structures (3) and (4) are actually not bound by 
a factor of 3 and 10, respectively, or that we underestimate 
the total gravitating mass on scales larger than 30 Mpc (cf. 
Fig. 9). 

We have used our cluster mass estimates (within -R500) 
to assess the probability that the various structures within 
the Shapley Supercluster arose from gaussian fluctuations 
in the primordial power spectrum. Adopting the CDM-likc 
power spectrum of Peacock & Dodds (1994) and fio = 1, 
we find that the supercluster core, which has an overdensity 
of between 4 and 11, is a fluctuation exceeding 3a which is 
approaching the turnaround point. If fio = 0.3, which al- 
lows the baryon fraction in the clusters to equal that of the 
rest of the Universe, then the significance of the superclus- 
ter core approaches 3<r with respect to the power expected 
on the same scale. These significance levels increase to 3.7 
and 3.5 a (for fio = 1 and 0.3, respectively) if the mass es- 
timates are extrapolated to i?200- If the significance of the 
fluctuations are required to be < 3<r, then fib > 0.095. Oth- 
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erwise, non-gaussian fluctuations need to be considered, or 
the power spectrum normalization to be increased on super- 
cluster scales. 

Further data are needed to make this extrapolation with 
confidence, but the probabilities implied for the structures 
emphasize the remarkable nature of, at least, the core of the 
Shapley Supercluster. 



ACKNOWLEDGEMENTS 

We are grateful to Harald Ebeling for his IDL support and 
to the members of the IoA X-ray Group. It is a pleasure to 
thank our referee, Gianni Zamorani, who has improved this 
work with his suggestions. ACF acknowledges the support of 
the Royal Society and DAW that of PPARC. This research 
has made use of data obtained through the High Energy As- 
trophysics Science Archive Research Center Online Service, 
provided by the NASA-Goddard Space Flight Centre. 



1989, ApJS, 70, 1 (ACO) 
: Szalay, A.S., 1986, ApJ, 



REFERENCES 

Abell, G., Corwin, H.G. & Olowin, R., 
Bardeen, J.M., Bond, J.R., Kaiser, N. 
304, 15 

Bardelli, S., Zucca, E., Vettolani, G., Zamorani, G., Scaramella, 
R., Collins, C. & MacGillivray, H.T., 1994, MNRAS, 267, 665 
Bardelli, S., Zucca, E., Malizia, A., Zamorani, G., Scaramella, R. 

& Vettolani, G., 1995, A&A, 305, 435 
Bartelmann, M. & Steinmetz, M., 1996, MNRAS, 283, 431 
Beers, T.C., Geller, M. & Huchra, J. P., 1982, ApJ, 257, 23 
Breen, J., Raychaudhury, S., Forman, W. & Jones, C, 1994, ApJ, 
424, 59 

Binney, J. & Tremaine, S., 1987, Galactic Dynamics, Princeton 

University Press 
Briel, U.G., Henry, J. P. & Bohringer, H., 1992, AA, 259, L31 
Carroll, S.M., Press, W.H. & Turner, E.L., 1992, ARAA, 30, 499 
Copi, C.J., Schramm, D.N. & Turner, M.S., 1995, Science, 267, 

192 

David, L.P., Slyz, A., Jones, C, Forman, W. & Vrtilek, S.D., 1993, 
ApJ, 412, 479 

David, L.P., Jones, C. & Forman, W., 1995, ApJ, 445, 578 
Day, C.S.R., Fabian, A.C., Edge, A.C. & Raychaudhury, S., 1991, 

MNRAS, 252, 394 
Dressier, A. & Shectman, S.A., 1988, AJ, 95, 985 
Edge, A.C, Stewart, G.C., Fabian, A.C. & Arnaud, K.A., 1990, 

MNRAS, 245, 559 
Einasto, M., Einasto, J., Tago, E., Dalton, G.B. & Andernach, 

H., 1994, MNRAS, 269, 301 
Evrard, A.E., Metzler, C.A. & Navarro, J.F., 1996, ApJ, 469, 494 
Fabian, A.C, Hu, E.M., Cowie, L.L. & Grindlay, J., 1981, ApJ, 

248, 47 
Fabian, A.C 
Gunn, K.F. 
Hata, N Stcig-m an. G 

ph/S 603087 



1991, MNRAS, 253, 29p 
Thomas, P.A., 1996, MNRAS, 281, 1133 

Bludman, S. & Langacker, P., 1996, lastro- 



Henry, J. P., Briel, U.G. & Nulsen, P.E.J. , 1993, A&A, 271, 413 
Kaiser, N. & Davis, M., 1985, ApJ, 297, 365 
Klemola, A.R., 1969, AJ, 74, 804 
Lacey, C. & Cole, S., 1993, MNRAS, 262, 627 
Makino, N. & Suto, Y., 1993, PASJ, 45, L13 
Melnick, J. & Moles, M., 1987, Rev Mex AA, 14, 72 
Mo, H.J. & White, S.D.M., 1996, MNRAS, 282, 347 
Padmanabhan, T., 1993, Structure formation in the universe, 
Cambridge University Press 



Peacock, J. A. & Dodds, S.J., 1994, MNRAS, 267, 1020 
Peebles, P.J.E., 1980, The Large-Scale Structure of the Universe, 

Princeton University Press 
Peebles, P.J.E., 1993, Principles of Physical Cosmology, Prince- 
ton University Press 
Postman, M., Huchra, J. P. & Geller, M.J., 1992, ApJ, 384, 404 
Quintana, H., Ramirez, A., Melnick, J., Raychaudhury, S. & 

Slezak, E., 1995, AJ, 110, 463 
Raychaudhury, S., 1989, Nature, 342, 251 

Raychaudhury, S., Fabian, A.C, Edge, A.C, Jones, C. & Forman, 

W., 1991, MNRAS, 248, 101 
Scaramella, R., Baiesi-Pillastrini, C, Chincarini, C, Vettolani, 

G. & Zamorani, C, 1989, Nature, 338, 562 
Schechter, P., 1976, ApJ, 203, 297 
Schindler, S., 1996, AA, 305, 756 

Schramm, D.N. & Turner, M.S., 1996, Nature, 381, 193 

Shapley, H., 1930, Harvard Obs. Bull., 874, 9 

Snowden, S.L., McCammon, D., Burrows, D.N. & Mendenhall, 

J. A., 1994, ApJ, 424, 714 
Stark A. A. et al , 1992, ApJS, 79, 77 

Stewart, G.C, Fabian, A.C, Jones, C. & Forman, W., 1984, ApJ, 
285, 1 

Tytler, D., Fan, X.M. & Buries, S., 1996, Nature, 381, 207 
Vettolani, C, Chincarini, C, Scaramella, R. & Zamorani, G., 

1990, AJ, 99, 1709 
Voges. W. et al , 1996, AA, in press 

Walker, T.P., Steigman, C, Schramm, D.N., Olive, K.A. & Kang, 

H. S., 1991, ApJ, 376, 51 

White, D.A. & Fabian, A.C, 1995, MNRAS, 273, 72 

White, D.A., Jones, C. & Forman, W., 1996, MNRAS, submitted 

White, S.D.M. & Frenk, C.S., 1991, ApJ, 379, 52 

White, S.D.M., Navarro, J.F., Evrard, A.E. & Frenk, C.S., 1993, 

Nature, 366, 429 
Zabludoff, A.I., Huchra, J. P. & Geller, M.J., 1990, ApJS, 74, 1 
Zucca, E., Zamorani, C, Scaramella, R. & Vettolani, C, 1993, 

ApJ, 407, 470 



APPENDIX A: A RECIPE TO EXTRAPOLATE 
THE LINEAR OVERDENSITY WHEN fi < 1 

In order to calculate the linear overdensity So at the present 
time and compare it with the observed non-linear overden- 
sity 8 both in an open Friedman-Robertson- Walker universe 
and in an Einstein-de Sitter universe with cosmological con- 
stant equal to and adopting a spherical collapse model, we 
refer to sections 11 and 19 of Peebles (1980, hereafter Pe), 
section 8 of Padmanabhan (1993, Pa) and appendix A in 
Lacey & Cole (1993, LC). 

The collapsing region is described in the evolution of 
its physical radius r during the time t by the parametric 
equations 



r = A c (l - cos9) 
t-T c = B c (6-sin8) 
Al = GMBl 



(Al) 
(A2) 
(A3) 



where T c will be ignored for a pure growing-mode perturba- 
tion (cf. Pa eq. 8.21). 

In the case of Oo = 1, the combination of the 3 equa- 
tions (Al-3) with the exact solution of the evolution of the 
background universe (cf. Pa eq. 8.23) and the extrapolated 
overdensity So in the linear regime (i.e. small t) provides the 
evolution of a spherical overdense region (cf. Pa 8.32-34): 
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1 + 5 = 



9(6>-sin6>) 2 
2(1 - cos(9) 3 



*~ 2H ( 3°) 



-3/2 



(A4) 



(A5) 



'W^^ 1 - -^^- (A6) 

Of these quantities, the observable at present time to arc 
the non-linear overdensity 5 and the comoving radius R = 
r(t)(l + z) allowing the detemination analytically of Ro and 
numerically of 9 and So. Their values are quoted in Table 4. 

The case for fio < 1 is discussed by Pe and a useful sum- 
mary is provided by LC. We present here the corresponding 
equations to (A4-6). In a open universe, the background 
with respect to a collapsing region that evolves according to 
(Al-3) is described by the parametric equations (cf. Pe eq. 
19.12) 



rb = At (cosh 77 — 1) 
t = Bj, (sinh 77 — rj) 
At = GUBl. 



(A7) 
(A8) 
(A9) 



Applying the two conditions that we are comparing to a per- 
turbated overdense region and a section of the background 
universe containing the same mass M and following the same 
time t, we obtain the relations 



Al 



B c (6 - sm6) = B h (sinh 77 - 77). 



(A10) 
(All) 



Then, following LC for extrapolation under linear condi- 
tions, we can write the linear overdensity S as 



So = 



3 sinh J7(sinhj7 — 77) 
(cosh 77 — l) 2 



1 + 



sinh 77 — 77 



2/3 



So = -D t 



Y sinh 77 — 77 J 



2/3' 



(A12) 



Thus, a growing density contrast can be described using 
(A10-11) by 

r]_ _ (cosh 77- l) 3 (g-sin6>) 2 



1 + 6 = 



1 



Ho 



(sinh 77 — t?) 2 (1 — cos 6>) 3 

3/2 



(A13) 



2H (1-fio) 3 / 2 \3D, 



2 So 



'to 



o 7T ^ (0-smO) (A14) 



r(t) = Ro 



2 So 

3 A 



1 — cos ( 



Ro 



2(fi- 1 -l) (1 + 5)V3' 



(A15) 



which represent the equations (A4-6) for Q < 1. Further- 
more, it can be shown that the approximation fig —t 1, 
i.e. 770 — > 0, provided by Taylor expansion of the hyper- 
bolic functions in (A13), gives the factor 9/2 as required for 
the case of a fiat universe, i.e. (A4). To estimate the values 
quoted in Table 4, we calculate numerically 60 at the present 
time to, given S and 770 = cosh _1 (2f2 ( 7 1 — 1). 

Finally, the time of maximum expansion (i.e. the 
turnaround point) t mc can be estimated using (A5) and 
(A14), for ilo = 1 and less than 1 respectively, with 6 = it. 
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